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The  formation  of  a  SEI  layer  and  its  growth  cause  internal  resistance  increase  and  capacity  loss,  leading  to 
performance  degradation  of  lithium-ion  batteries.  In  order  to  comprehensively  investigate  the  effects  of 
SEI  growth  on  battery  performance,  a  one-dimensional  thermal-electrochemical  model  was  developed. 
This  model  is  equipped  with  a  growth  mechanism  of  the  SEI  layer  coupled  with  thermal  evolution,  based 
on  the  diffusional  process  of  the  solvent  through  the  SEI  layer  and  the  kinetic  process  at  the  interface 
between  the  solid  and  liquid  phases.  The  model  is  able  to  reveal  the  effects  of  diffusivity,  reaction  kinetics 
and  temperature  on  SEI  layer  growth  and  cell  capacity  fade.  We  show  that  depending  on  the  SEI 
thickness,  the  growth  can  be  kinetics-limited  or  diffusion-limited.  With  the  layer  becoming  thicker,  its 
growth  rate  slows  down  gradually  due  to  increased  diffusion  resistance.  The  SEI  layer  grows  faster  during 
charge  than  discharge  due  to  the  difference  in  the  electron  flux  through  the  SEI  layer  and  the  temper¬ 
ature  change  during  cycling.  Temperature  rise  due  to  reaction  and  joule  heating  accelerates  the  SEI  layer 
growth,  leading  to  more  capacity  loss.  Our  model  can  provide  insights  on  position-dependent  SEI  growth 
rate  and  be  used  to  guide  the  strategic  monitoring  location. 
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1.  Introduction 

Aging  is  a  significant  concern  in  Lithium-ion  Batteries  (LIBs), 
especially  in  the  application  to  electric  vehicles  (EVs)  and  hybrid 
EVs  (HEVs)  because  of  the  requirement  of  cell  durability  [1,2]. 
Battery  aging  is  caused  by  the  material  degradation  of  cell  com¬ 
ponents  and  the  introduction  of  new  components  by  side  reactions. 
It  is  well  known  that  electrolyte  decomposition,  which  is  linked  to 
the  side  reactions  at  electrode  surfaces,  yields  a  passivation  film 
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[3-5  .  This  film,  referred  as  the  solid  electrolyte  interface  (SEI), 
suppresses  further  lithium  corrosion  and  solvent  decomposition 
[6,7  .  Long-term  stability  of  the  SEI  would  have  a  prime  impact  on 
capacity  retention,  internal  resistance  and  cyclability  of  LIBs.  It  has 
been  reported  that  the  long-term  degradation  of  large-scale  LIBs  is 
caused  by  the  formation  and  continuous  growth  of  the  SEI  layer  8]. 
In  spite  of  numerous  research  efforts,  fundamental  understanding 
of  the  SEI  composition,  stability,  and  its  influence  on  battery  per¬ 
formance  is  still  controversial,  and  has  been  a  focus  of  many  re¬ 
searches  over  the  past  decade  [9-13]. 

Developing  appropriate  laboratory  aging  tests  to  analyze  the 
relation  between  SEI  formation  and  cell  aging  is  challenging.  SEI 
growth  involves  a  number  of  simultaneous  physicochemical  pro¬ 
cesses  and  their  interactions  [7,14-16  .  Theoretical  modeling  and 
simulations  are  very  useful  to  provide  insight  into  the  effects  of 
SEI.  Mathematical  models  have  been  developed  to  simulate  the 
growth  of  SEI  and  lithium  transport  in  the  compact  layer  of  SEI 
[9,17-19]. 

Advances  in  the  understanding  of  SEI  layer  formation  and 
growth  call  for  a  model  to  account  for  new  observations.  Research 
works  suggest  that  the  reductive  electrolyte  decomposition  is 
accompanied  by  irreversible  lithium  ion  consumption  at  the  elec¬ 
trode/electrolyte  interface.  The  products  from  decomposition  and 
side  reactions  build  up  and  form  the  SEI  layer  [7  .  Some  recent 
studies  show  that  the  SEI  layer  on  the  negative  electrode  material  is 
composed  of  two  distinct  layers:  a  compact  polycrystalline  or 
heteromicrophase  layer  of  inorganic  products  (e.g.  Li2C03,  LiF,  Li20, 
etc)  and  a  porous  amorphous  outer  layer  composed  of  partially 
reduced  organic  products  (e.g.  (CH20C02Li)2,  ROLi,  R0C02Li,  etc) 
[20-26].  It  is  suggested  that  the  outer  layer  of  the  SEI  layer,  whose 
pores  are  filled  with  electrolyte,  controls  the  progress  of  SEI  growth, 
and  consequently  the  extent  of  capacity  fade  and  internal  resis¬ 
tance  rise.  The  outer  layer  can  be  a  couple  of  order  of  magnitudes 
thicker  than  the  inner  layer,  while  at  the  same  time  the  outer  layer 
is  more  vulnerable  and  unstable  [23,27-29  .  A  recent  work  inves¬ 
tigated  the  diffusion  phenomena  in  the  SEI  porous  layer  and  its 
growth  [13  .  However,  the  model  was  for  a  single  particle. 


SEI  growth  depends  on  the  electron  current  flow  and  lithium 
diffusion,  which  vary  at  different  locations  in  the  electrode.  A  cell- 
level  physics-based  SEI  growth  model  is  necessary  to  reveal  this 
spatial  distribution  of  SEI  and  evaluate  its  effect  on  cell  performance. 
Moreover,  it  is  important  to  integrate  the  temperature  effect  in  the 
cell-level  model  since  it  directly  influences  the  material  properties 
such  as  diffusivity  and  side  reaction  rates.  Temperature  is  also  of 
interest  for  the  study  of  thermal  runaway  which  is  closely  related  to 
the  safety  of  LIBs.  Several  cell-level  thermal-electrochemical  models 
were  developed  to  investigate  the  temperature  distribution  and  its 
effect  on  cell  performance,  however,  SEI  layer  was  either  neglected 
or  simplified  as  constant  resistance  in  the  battery  system  [30-32]. 
Thermal  effect  majorly  impacts  the  side  reaction  for  SEI  formation 
and  growth  in  the  cell,  while  the  thermal  behavior  of  LIBs  is  in  turn 
determined  by  the  chemical  and  electrochemical  processes  occur¬ 
ring  inside  the  cell  during  cycling  [30  .  Therefore,  the  SEI  growth 
analysis  needs  to  be  coupled  with  the  thermal  effect. 

In  this  paper  a  cell-level  physics-based  model  is  developed  with 
the  consideration  of  multiple  key  electrochemical  processes, 
particularly  the  side  reaction  for  continuous  SEI  growth  and  the 
thermal  effect.  The  model  has  the  potential  to  provide  insights  into 
the  mechanisms  related  to  capacity  and  power  loss,  as  well  as  to 
reveal  their  effects  on  electric  potential  and  Li  ion  concentration 
distribution.  By  considering  solvent  transport  through  the  SEI 
porous  layer  [13  ,  the  model  can  capture  lithium-ion  diffusion  and 
side  reaction  kinetics  inside  the  layer.  The  model  would  also  enable 
us  to  predict  position-dependent  SEI  growth,  providing  important 
guidance  for  monitoring  locations.  In  addition,  this  model  considers 
the  effect  of  SEI  compact  layer  on  lithium  ion  diffusion.  SEI  layer 
formation  and  growth  are  quite  temperature-dependent.  By  incor¬ 
poration  of  thermal  energy  conservation,  the  interaction  between 
SEI  layer  growth  and  temperature,  as  well  as  their  coupled  effects  on 
battery  performance  can  be  described.  The  model  also  considers 
double-layer  capacity  [33,34  ,  which  allows  future  electrochemical 
impendence  spectroscopy  (EIS)  simulation.  EIS  simulation  is  a 
powerful  tool  to  interpret  Li  ion  and  electron  migration  mecha¬ 
nisms,  interfacial  phenomenon,  and  failure  mechanisms  [35  . 
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Fig.  1.  A  schematic  showing  the  SEI  layer,  double-layer  capacity  and  thermal  effect  in  a  lithium  ion  battery. 
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2.  Model  development 

Fig.  1  shows  a  schematic  of  the  SEI  layer,  double-layer  capacity 
and  thermal  effect  in  a  typical  cell  consisting  of  the  current  col¬ 
lector,  anode,  cathode,  separator,  and  electrolyte  solution.  The  SEI 
layer  on  the  anode  is  much  thicker  compared  to  that  on  the  cath¬ 
ode,  which  causes  more  rise  in  cell  internal  resistance  and  more 
capacity  loss  [36].  Therefore,  here  we  focus  on  the  SEI  layer  growth 
on  the  anode.  SEI  formation  and  growth  reaction  occurs  continu¬ 
ously  between  the  lithiated  graphite  and  the  electrolyte  solvent.  In 
the  meanwhile,  the  growing  SEI  layer  gradually  consumes  lithium 
and  causes  cell  impedance  increase  during  cycling.  In  this  study  the 
electrolyte  solvent  mixture  contains  ethylene  carbonate  (EC).  The 
SEI  is  considered  to  be  made  of  (CH20C02Li)2,  which  has  been 
shown  to  be  a  major  component  of  the  SEI  porous  layer  [20-26]. 
The  SEI  formation  process  can  be  described  by  the  following  re¬ 
actions  [13  . 

EC  +  e~  (graphite)  ->  EC~  (1) 

EC-  +EC-->(02C0)-(CH2)2(02C0)-  +C2H4T  (2) 


the  SEI  thickness  at  time  t  and  ctsei  is  the  SEI  ionic  conductivity.  The 
cathode  SEI  resistance  is  assumed  zero  since  it  is  much  thinner.  In 
Eq.  (6)  the  first  term  accounts  for  the  current  flux  due  to  conduc¬ 
tivity,  and  the  second  term  accounts  for  the  effect  of  potential 
variation  due  to  concentration  variation  of  the  salt  in  the  electro¬ 
lyte,  where  cr|ff  is  the  effective  conductivity  of  the  electrolyte,  F  is 
the  Faraday's  constant,  R  is  the  universal  gas  constant,  T  is  tem¬ 
perature,^  is  the  activity  coefficient  of  salt  in  the  electrolyte  phase, 
c2  is  the  lithium  concentration  in  the  electrolyte,  and  is  the 
transference  number  of  lithium  ion  with  respect  to  the  solvent. 
Equation  (7)  gives  the  intercalation  current  density  calculated  from 
the  solid  phase,  where  is  the  effective  conductivity  of  the 
electrode  and  the  current  density  in  the  solid  phase  is 
i\  =  —  (00-[  /dx). 

The  material  balance  for  the  electrolyte  in  the  solution  phase  is 
given  by 

9(T2c2)  _  9  f neffdc2\  ,  asi\ oc(A  i.o\ 

~~dF~  ~  dx  V  2 

as^  (a  i-0\  9(01  -  02  -  Ef  -  flsElhoc)  fQ\ 

+  Tl1“t+J - at -  (8) 


EC-  +  £C  +  e-(graphite)^(02C0)“(CH2)2(02C0)“  +  C2H4T 

(3) 

(02C0)-(CH2)2(02C0)- 

+  2Li+— >Li+(02C0)-(CH2)2(02C0)-Li+  (4) 

The  double-layer  capacity  represents  an  additional  mechanism 
to  account  for  energy  storage  and  the  response  to  rapid  current 
pulses.  A  modification  to  the  regular  current  balance  equations  is 
needed  to  consider  the  charging  of  double  layers  in  both  anodes 
and  cathodes  [33,34].  The  modification  also  appears  in  the  mass 
balance  equation  for  the  electrolyte,  since  it  includes  the  diver¬ 
gence  of  the  current  as  a  homogeneous  source  term.  The  thermal 
effect  will  be  considered  by  a  microscopic-macroscopic  modeling 
approach. 


where  r2  is  the  electrode  porosity  or  electrolyte  phase  volume 
fraction,  and  D|ff  is  the  effective  Li  diffusion  coefficient  in  the 
electrolyte.  The  last  term  on  the  right  of  Eq.  (8)  accounts  for  the 
extra  salt  storage  variation  over  time  in  the  solution  phase  due  to 
the  double  layer  capacity  effect. 

The  diffusion  of  Li  ions  in  an  active  solid  particle  is  governed  by 
the  Fields  law 


9Ci  n  f 9 2C\  2  OcA 

at  1  ^  a r2  +  r  ar  J 


where  C\  is  the  concentration  of  Li  in  the  particle,  r  is  the  radial 
distance  in  the  spherical  coordinate  and  Di  is  the  diffusion  coeffi¬ 
cient  of  Li  in  the  solid. 

The  solid  and  solution  phases  are  related  by  the  boundary 
condition, 


2.1.  ID  cell  model 


The  charge  conservation  in  the  electrode  and  electrolyte  phase 
can  be  described  by 
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Equation  (5)  gives  the  intercalation  current  density  in  the 
porous  electrode,  where  as  is  the  active  surface  area  per  unit  elec¬ 
trode  volume  and  ijoc  is  the  local  current  density.  The  first  them  on 
the  right  of  Eq.  (5)  accounts  for  the  divergence  of  the  current  flux  in 
the  electrolyte,  where  i2  is  the  current  density  in  the  electrolyte 
phase.  The  second  term  accounts  for  the  extra  current  due  to  the 
double  layer  capacity  effect,  where  C  is  the  double-layer  capaci¬ 
tance,  0i  is  the  electric  potential  of  the  solid  phase,  02  is  the  electric 
potential  of  the  liquid  phase,  U  is  the  open-circuit  potential  that  is 
dependent  on  the  state  of  change  (SOC),  and  Rsei  is  the  SEI  layer 
resistance.  In  this  study  we  focus  on  the  SEI  layer  on  the  anode.  The 
anode  SEI  resistance  is  estimated  by  RseM  =  W/asEi,  where  <5(t)  is 


(10) 

r=Rs 

where  in  is  the  current  density  of  Li  intercalation  into  or  dein¬ 
tercalation  from  graphite  and  Rs  is  the  radius  of  the  active  material 
particle. 

The  intercalation  current  at  cathode  can  be  described  by  the 
Bulter-Volmer  equation, 


in  =  -FD- 
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(11) 


where  ko  is  the  rate  coefficient  for  lithium  intercalation  reaction  at 
cathode,  (3  is  the  charge-transfer  coefficient,  CitS  is  the  Li  concen¬ 
tration  on  the  surface  of  the  active  material  particle,  and  ci>max  is 
the  maximum  possible  Li  concentration  in  the  active  material. 


2.2.  SEI  growth  modulus 

Here  we  extend  SEI  growth  on  the  surface  of  a  single  particle  to  a 
thermal-electrochemical  ID  cell-level  model.  The  electrolyte 
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solution  is  LiPF6  in  EC/DMC.  The  compact  SEI  layer  is  assumed  to  be 
1  nm  thick  (Lsei(O)  =  1  nm),  which  also  accounts  for  an  initial  SEI 
layer  resistance  of  Ksei(O)  =  1Q  cm2.  The  porous  layer  is  assumed  to 
grow  above  the  compact  layer.  The  dilute  solution  theory  can  be 
used  to  describe  the  flux  of  charged  species  in  the  SEI  layer  [37]. 
Using  the  Bulter-Volmer  equation  to  describe  the  lithium  inter¬ 
calation/deintercalation  current  in  the  active  materials  occurring  at 
the  interface  between  the  graphite  and  the  SEI  layer,  we  have 


'Li  =  W<0C',  s(Cl  max  -  Cl  s)’  ffc\  T 


exp 


(1  -  0)F 

RT 


(01  —  02  —  ^ 


-  KsElhoc) 


-  exp 


1 6F 


-  (01  -  02  —  U  -  ^SElhoc) 


(12) 


where  /<o  is  the  rate  coefficient  for  lithium  intercalation  reaction  at 
anode  and  (3  is  the  charge-transfer  coefficient. 

The  rate  determining  step  is  considered  to  be  the  radical  anion 
formation  shown  in  Eq.  (1).  The  side  reaction  takes  place  at  the 
same  interface  of  the  lithium  intercalation/deintercalation  reaction, 
and  the  kinetic  expression  for  this  side  reaction  is  in  the  form  of  a 
Tafel-like  expression, 


(13) 

where  /<sei  is  the  rate  constant  of  side  reaction  for  SEI  formation,  cec 
is  the  solvent  concentration  in  the  SEI  layer,  and  /3Sei  is  the  charge- 
transfer  coefficient.  Considering  simultaneous  lithium  intercala¬ 
tion/deintercalation  reaction  and  side  reaction,  we  can  write  the 
local  current  density  at  anode, 

hoc  =  *Li  +  *side  (14) 

Now  we  look  at  the  solvent  transport  in  the  SEI  porous  layer.  The 
rate  of  the  electrolyte  solution  concentration  change  inside  the  SEI 
porous  layer  is  compensated  by  the  gradient  of  electrolyte  solution 
flux  density,  Nec,  which  can  be  described  by 


hide  =  —FksEl cEc\  exP 


@se\F 

RT 


(01  -  02  -  El  hoc) 


dggc  _  _9Nec 
at  dr 


(15) 


Considering  both  diffusion  and  convection,  the  flux  density  is 
given  by 


Nec  =  -Dec-^t  +  vCec  (16) 

where  Dec  is  the  diffusion  coefficient  and  v  is  the  velocity  of  the 
solvent  relative  to  the  SEI.  The  growth  rate  of  the  SEI  layer  due  to 
the  accumulation  of  lithium  ethylene  dicarbonate  is  related  to  the 
side  reaction  rate  by 


d<5  _  *side  MsEI 
dt  2 F  pSEI 


(17) 


where  Msei  is  the  molecular  weight  of  the  SEI  and  psei  is  the  density 
of  the  SEI. 

Substitution  of  Eq.  (16)  to  Eq.  (15)  gives  the  material  balance  in 
the  SEI  porous  layer 
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At  the  boundary  between  the  SEI  and  the  particle  surface,  the 
solvent  flux  matches  the  side  reaction  current,  which  gives 


NEc\r=Rs  = 


*side 
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(19) 


At  the  interface  between  the  SEI  and  the  solvent,  the  concen¬ 
tration  is  given  by 


cEc\r=Rs+5  =  £SEIcEC,0  (20) 

where  rsEi  is  the  SEI  porosity  and  cec, o  is  the  solvent  concentration 
in  the  bulk  of  the  electrolyte. 

The  boundary  conditions  for  the  cell  level  model  are 


dC-[ 

—  0 

0C! 

hi 

dC2 

dC2 

dr 

r=0 

dr 

r=Rs 

FDh 

dX 

x=0 

~  dX 

(21) 


Fig.  2.  Coupling  of  the  models. 


Table  1 

Values  of  parameters  used  in  the  simulations. 


Parameters 

Anode 

Cathode 

Separator 

SEI 

Active  surface  area  per  unit  electrode  volume  (m-1),  as 

113040 

111375 

Cell  area  (m2),  A 

2.4  x  10-3 

Maximum  solid  phase  Li  concentration  (mol  m-3),  climax 

26390 

22860 

Initial  electrolyte  salt  concentration  (mol  m  3),  c2(0) 

2000 

EC  solution  phase  concentration  (mol  m-3),  Cce, o 

4541 

Double-layer  capacitance  (F  m-2),  Cc 

0.1 

0.1 

Heat  capacity  (J  kg  1  K-1),  Cp 

700 

700 

700 

Solid  phase  Li  diffusion  coefficient  (m2  s-1),  Dia 

3.9  x  10-14 

1  x  10  13 

Electrolyte  phase  Li  diffusion  coefficient  (m2  s-1),  D2 

7.5  x  10-11 

7.5  x  10-11 

7.5  x  10-11 

Effective  electrolyte  phase  Li  diffusion  coefficient  (m2  s-1), 

Df  =  D2(c2f 

EC  diffusion  coefficient  (m2  s-1),  DEC 

1  x  10 

-13 

Activation  energy  for  Li  diffusion  in  solid  phase  (kj  mol-1),£act,Di 

4 

20 

Activation  energy  for  Li  diffusion  in  electrolyte  phase  (kj  mol_1),Eact,D2 

10 

Activation  energy  for  exchange  current  density  (kj  mol-1),£actj/co 

30 

30 

Activation  energy  for  SEI  reaction  rate  (kj  mol _1),£act,/<sEib 

200 

Activation  energy  for  electrolyte  conductivity  (kj  mol-1),  £ac t,ff2 

20 

Electrolyte  activity  coefficient, /+a 

1 

Faraday's  constant  (C  mol-1),  f 

96485.3 

Heat  transfer  coefficient  (W  m-2  K-1),  h 

5 

5 

Reaction  rate  coefficient  (m4-3*3  s-1  mol*3-1),  ko 

2  x  10  6 

2  x  10-6 

SEI  reaction  rate  coefficient  (m  s-1),  kSEi 

1.37  x 

10-5 

Electrode  thickness  (pm),  L 

100 

183 

52 

Initial  SEI  thickness  (nm),  Lsei(O) 

1 

SEI  molar  mass  (kg  mol-1),  Ms ei 

0.162 

Bruggeman  porosity  exponent,  p 

1.5 

1.5 

Gas  constant  (J  mol-1  K-1),  R 

8.314 

Particle  radius  (pm),  Rs 

12.5 

8.5 

SEI  initial  resistance  (Q  cm2),  Rsei(O) 

1 

Initial  state  of  charge,  SOC 

0.56347 

0.1706 

Li  transference  number, t9 

0.363 

0.363 

0.363 

Environmental  temperature  (K),  T0 

298.15 

Open  circuit  potential  (V),  U 

U 

U-  =  -0.16  +  1.32  exp(3.  Ox)  +  10.0  exp( -2000.  Ox)  forLixC6,  ^  =  C\s/c\mAX  atanode 

U+  =4.19829  +  0.0565661  tanh[-14.5546y  +  8.60942]  -  0.0275479[l/(0.998432  _y)°  492465 

-  1.90111]  -  0.157123  exp (0. 0473 8y 

8)  +  0.810239  exp[-40(y  - 

0.133875)]  for  LiyMn204, 

y  =  Ci,s/ci, max  at  cathode 

Entropy  change  of  the  electrode  reaction  (V  K-1),  dU/dT1 

0 

0 

Charge  transfer  coefficient,  (3 

0.5 

0.5 

SEI  charge  transfer  coefficient,  /3Sei 

0.5 

Solid  phase  volume  fraction, £+ 

0.471 

0.297 

Electrolyte  phase  volume  fraction,  e2 

0.357 

0.444 

1 

SEI  Porosity,  £SEif 

0.05 

Thermal  conductivity  (W  m-1  K-1),  A 

5 

5 

1 

Mass  density  (kg  m-3),  p 

2.5  x  106 

1.5  x  106 

1.2  x  106 

SEI  density  (kg  m-3),  pSEi 

1.69  x 

103 

Solid  phase  conductivity  (S  m-1),^3 

100 

3.8 

Effective  solid  phase  conductivity  (S  m-1),  cr^ff 

<+  =  <7l(£l)p 

Electrolyte  phase  ionic  conductivity  (S  m-1),  a2 

<7 2  =  1.0793  X  10-2  +6.7461  x  10-4c2  -  5.2245  x  10-7(c2)2  +  1.3605  x  10-1°(c2)3  -  1.1724 

X  10-+c2)4 

Effective  electrolyte  phase  ionic  conductivity  (S  m-1),  cr|ff 

=  G2(£2)P 

Effective  electrolyte  phase  diffusional  conductivity  (A  m-1),  o 

eff  _  2  RTaf  n  fo^  1  d  In  f+\ 

a2D  ~  F  V  1  C+  H  1  +  d  In  c2  1 

SEI  ionic  conductivity  (S  m-1),  aSE\ 

5  x  10 

6 

a  Data  from  Ref.  [41  ]. 
b  Data  from  Ref.  [12]. 
c  Data  from  Ref.  [34]. 
d  Data  from  Ref.  [39]. 
e  Data  from  Ref.  [17]. 
f  Data  from  Ref.  [13]. 
s  Assumed. 


486  L.  Liu  et  al.  /  Journal  of  Power  Sources  268  (2014)  482-490 


L.  Liu  et  al.  /  Journal  of  Power  Sources  268  (2014)  482-490 


487 


eff001 


9x 


eff90i 

I  9  0i 

d</>i 

-  ffi  dx 
x=0  ox 

x=L  A?  dX 

x=L _ 

x — L_  -\ -L 


sep 


0. 


902 


9x 


902 


x=0 


9x 


0 


x=L 


(22) 


where  L  is  the  thickness  of  the  cell,  /  is  the  cell  current,  A  is  the 
cross-section  area  of  the  cell,  L_  is  the  thickness  of  the  anode,  and 
Lsep  is  the  thickness  of  the  separator. 


2.3.  Thermal  modulus 


U) 


Time  (s)  XI 03 


The  thermal  energy  conservation  can  be  expressed  by 


d(pCpT) 
9 1 


VAVT  +  q 


(23) 


where  p  is  the  mass  density,  Cp  is  the  heat  capacity,  and  A  is  the 
thermal  conductivity  of  the  cell  component.  The  heat  generation 
rate  per  unit  volume,  q ,  can  be  expressed  by  Refs.  [38,39], 

q  =  asi\oc  ^0i  —02  —  U  —  ftsEihoc  +  T  +  <rfff(V0i)2 

+  (V02)2  +  ^2DVln  C2  •  V02  (24) 

where  ct^d  is  th e  effective  solution  phase  diffusional  conductivity. 
Here,  the  first  term  of  the  right-hand  side  represents  the  heat 
sources  due  to  charge  transfer  at  the  electrode/electrolyte  in¬ 
terfaces  and  the  entropy  change.  The  last  two  terms  arise  from  the 
joule  heating  in  the  solid  active  and  electrolyte  phases,  respectively. 
The  boundary  conditions  are  given  at  x  =  0  and  x  =  L 
by  -AVT  =  /i(T-To),  where  To  is  the  environmental  temperature, 
and  h  is  the  heat  transfer  coefficient. 

The  temperature-dependent  physicochemical  properties  are 
descripted  by  the  Arrhenius  equation, 


4>  =  4>ref  exp 


(25) 


where  <f>  is  a  general  variable  that  represents  the  diffusion  coeffi¬ 
cient  of  a  species,  the  conductivity  of  the  electrolyte,  or  the  ex¬ 
change  current  density  of  an  electrode  reaction,  etc.,  with  the 
subscript  ref  indicating  the  value  at  a  reference  temperature.  Eact  is 
the  corresponding  activation  energy. 

We  developed  our  multi-scale  and  multi-physics  model  using 
the  COMSOL  Multiphysics  4.2a  software.  Fig.  2  shows  the  coupling 
among  the  components.  The  ID  porous  electrode  model  is 
composed  of  the  anode  region  (0  <  x  <  L_),  the  separator  region 
(L_  <  x  <  L_  +  Lsep),  and  the  cathode  region  (L_  +  Lsep  <  x  <  L).  To 
couple  the  ID  model  with  the  Li+  diffusion  in  the  spherical  solid 
particle,  we  add  a  vertical  axis,  r,  to  form  a  pseudo  2D  model.  This 
vertical  axis  denotes  the  radial  direction  of  the  solid  particles.  In  the 
two  dimensional  (x,  r)  space,  the  Li+  diffusion  coefficient  in  the  x 
direction  is  set  to  several  orders  of  magnitude  lower  than  that  in 
the  r  direction  to  accomplish  the  radial  diffusion  descried  by  Eq.  (9). 
The  current  flux  calculated  from  the  ID  model  is  projected  to  the 
lower  r  =  Rs  boundary  in  Fig.  2,  to  solve  the  Li+  concentration,  C\,  in 
the  solid  particle  phase.  The  solved  concentration  is  then  projected 
back  to  the  1 D  model  to  provide  the  surface  concentration  of  the 
solid  particle  phase.  Similarly,  the  SEI  growth  is  performed  in  the 
pseudo  2D  configuration  and  then  projected  back  to  the  ID  model 
to  provide  the  resistance  of  the  SEI  layer. 


Fig.  3.  SEI  layer  growth  and  its  resistance  reach  a  stable  state  after  7  cycles  at  the 
charge/discharge  rate  of  1C.  SEI  layer  grows  faster  during  charge  than  during  discharge. 


3.  Results  and  discussion 

The  values  of  the  parameters  used  in  the  simulations  are  listed 
in  Table  1.  Our  model  has  the  transient  feature  and  can  provide  the 
capability  to  evaluate  SEI  growth  during  charge  and  discharge  cy¬ 
cles.  An  initial  SEI  layer  thickness  of  1  nm  was  assumed  to  account 
for  the  SEI  compact  layer.  The  compact  layer  was  assumed  to  form 
during  the  first  charge  of  the  cell  to  a  SOC  of  56%.  The  porous  SEI 
layer  formed  and  grew  outside  the  compact  layer.  The  growth  of  the 
porous  layer  is  very  sensitive  to  the  diffusion  process  and  the  side 
reaction  rate.  Fig.  3  shows  the  SEI  layer  growth  and  the  increase  in 
layer  resistance  during  the  first  seven  cycles  at  1C  charge/discharge 
rate.  After  the  initial  discharge  process,  the  SEI  layer  grows  by 
3.8  nm.  During  the  next  charge  and  discharge  cycles,  the  SEI  layer 
grows  by  2  nm  and  0.6  nm,  respectively.  The  growth  continuously 
slows  down  during  further  cycling.  After  seven  cycles  the  thickness 
of  the  SEI  porous  layer  increases  to  around  9  nm  while  its  resistance 
rises  to  18  Q  cm2.  Our  simulations  suggest  that  the  SEI  porous  layer 
grows  quickly  toward  the  equilibrium  thickness.  This  trend  agrees 
with  the  reports  9,13,17,40]  showing  that  SEI  growth  tends  to  slow 
down  due  to  possible  difficulty  in  electrolyte  solution  diffusion  or 
high  electrical  resistance  resulted  from  the  SEI  layer.  Fig.  3  also 
shows  that  the  SEI  layer  grows  faster  during  charge  than  during 
discharge,  which  is  consistent  with  experimental  observations  [9]. 
The  growth  rate  is  determined  by  the  electron  flux  through  the  SEI 
layer.  During  charge,  the  injection  of  electrons  into  the  anode 
makes  it  easier  to  leave  the  graphite  and  react  with  electrolyte  and 
lithium  ions  to  form  SEI.  The  difference  of  growth  rate  during 
charge  and  discharge  may  also  be  related  to  the  temperature.  Fig.  4 


Time  (s) 

Fig.  4.  A  higher  C-rate  leads  to  more  temperature  increase.  At  the  same  C-rate,  the 
temperature  during  charge  is  a  litter  higher  than  that  during  discharge. 
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Fig.  7.  Effect  of  environmental  temperature  on  SEI  growth. 
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shows  that  the  temperature  is  slightly  higher  during  charge  than 
discharge,  and  the  difference  becomes  more  noticeable  as  the  C- 
rate  becomes  larger.  SEI  grows  faster  at  a  higher  temperature. 

It  is  important  to  understand  the  relationship  among  diffusivity, 
side  reaction  rate,  and  SEI  layer  growth  characteristics.  When  the 
diffusion  coefficient  is  small,  the  system  becomes  diffusion-limited. 
Fig.  5  shows  the  SEI  layer  growth  with  different  diffusion  co¬ 
efficients  for  the  electrolyte  in  the  SEI  layer.  With  the  diffusion 
coefficients  increased  by  2  orders  of  magnitude,  which  is  consistent 
with  the  possible  range  reported  in  the  literatures  [13,41  ,  the  SEI 
layer  grows  5  times  as  thick.  With  a  lower  diffusion  coefficient,  the 
SEI  layer  tends  to  reach  its  stable  state  faster  but  at  a  thinner 
thickness.  For  example,  the  SEI  layer  stops  growth  at  cycle  2  and 
reaches  a  thickness  of  4  nm  with  Dec  =  0.1  Do,  where  Do  is  the 
baseline  value  given  in  able  1.  If  the  diffusion  coefficient  rises  by  2 
orders  of  magnitude,  the  SEI  layer  reaches  a  stable  state  at  around 
cycle  7  or  8  with  a  larger  thickness  of  20  nm.  With  a  higher  diffusion 
coefficient,  the  electrolyte  solution  penetrates  through  the  SEI 
porous  structure  more  easily  to  reach  the  reaction  sites.  Therefore, 
the  SEI  layer  grows  faster  and  thicker.  The  SEI  layer  growth  also 
depends  on  its  kinetic  reaction.  With  a  large  diffusion  coefficient, 
the  system  becomes  more  kinetics-limited.  Fig.  6  shows  the  effect 
of  side  reaction  rate  constant  on  SEI  growth.  /<o  is  the  baseline  value 
given  in  Table  1.  The  reaction  rate  can  vary  2-3  orders  of  magnitude 
depending  on  the  system  [9,13].  With  a  faster  reaction  rate,  the  SEI 
layer  grows  faster  and  thicker,  but  takes  a  longer  time  to  reach  the 
stable  thickness.  The  overall  trend  is  that  the  SEI  grows  quickly  in 
the  first  few  cycles.  The  initial  growth  of  SEI  is  kinetics-limited.  The 
growth  rate  gradually  slows  down  due  to  the  rising  diffusional 
resistance  from  the  layer  thickness,  making  the  system  shift  toward 
diffusion-limited. 


The  performance  of  LIBs  can  be  significantly  influenced  by 
temperature  and  the  thermal  behavior,  which  is  also  closely 
related  to  battery  safety.  While  the  thermal  process  can  accelerate 
the  electrochemical  and  chemical  processes  in  LIBs,  they  also  cause 
problems  such  as  more  internal  resistance  increase  and  capacity 
loss  due  to  faster  decomposition  of  the  electrolyte  and  SEI  growth. 
Fig.  7  shows  the  temperature  effect  on  SEI  growth.  With  a  higher 
environmental  temperature,  the  side  reaction  kinetics  and  diffu¬ 
sion  accelerate,  resulting  in  a  thicker  SEI  layer.  At  45  °C,  the  SEI 
thickness  increases  by  more  than  15%  in  the  first  8  cycles 
comparing  to  that  at  room  temperature.  Previous  studies  have 
suggested  that  a  higher  temperature  will  accelerate  the  diffusion 
process  [39,42  .  Here  we  show  quantitatively  how  temperature 


Fig.  8.  A  higher  environmental  temperature  leads  to  more  ohmic  loss  and  capacity  loss 
due  to  SEI  growth. 
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Fig.  10.  Effect  of  C-rate  on  SEI  growth. 


Fig.  11.  Discharge  profile  at  1C  with  and  without  SEI  growth. 


affects  the  SEI  growth.  The  resulting  increase  of  internal  resistance 
due  to  a  thicker  SEI  leads  to  more  ohmic  and  capacity  loss,  as 
shown  in  Fig.  8.  We  have  investigated  the  increase  of  cell  average 
temperature  as  a  result  of  different  boundary  conditions,  as  shown 
in  Fig.  9.  With  a  higher  thermal  transfer  coefficient,  the  increase  in 
cell  average  temperature  is  negligible.  However,  as  the  cooling 
condition  becomes  worse,  the  cell  temperature  rises  dramatically. 
As  we  have  shown,  higher  temperature  results  in  more  ohmic  loss 
and  capacity  loss.  Therefore,  temperature  control  is  critical  in  cell 
performance.  C-rate  is  another  important  parameter  in  deter¬ 
mining  SEI  layer  growth  and  cell  performance.  Fig.  10  shows  the 
effect  of  C-rate  on  SEI  growth.  The  SEI  layer  grows  almost  4  times 


XI  O'13 


Fig.  12.  Position-dependent  SEI  growth.  The  growth  rate  is  higher  at  the  anode/ 
separator  interface. 


as  thick  when  the  C-rate  is  increased  from  0.5C  to  1.5C,  which 
implies  that  the  SEI  layer  growth  is  highly  dependent  on  the  cur¬ 
rent  flow  density  inside  the  cell.  The  increase  of  internal  resistance 
and  capacity  loss  due  to  SEI  growth  will  result  in  the  degradation 
in  cell  performance.  Fig.  11  shows  the  comparison  of  discharge 
profile  at  1C  with  and  without  SEI  growth.  The  SEI  layer  resistance 
increases  to  around  60  Q  cm2  after  10  cycles,  which  causes  about 
4%  capacity  loss.  The  voltage  drop  mainly  comes  from  the 
increased  cell  internal  resistance  associated  with  the  SEI  growth. 
The  SEI  growth  depends  on  electron  current  flow  and  lithium 
diffusion,  which  varies  at  different  locations  in  the  electrode.  A 
cell-level  SEI  growth  model  can  reveal  position-dependent  SEI 
growth,  which  can  be  used  as  guidance  for  the  location  of  moni¬ 
toring.  Fig.  12  shows  the  SEI  growth  rate  along  the  anode  thickness 
direction.  The  SEI  layer  grows  faster  at  the  interface  between  the 
anode  and  the  separator  (x  =  L_)  than  at  the  current  collector  side 
(x  =  0).  The  non-uniform  SEI  growth  can  result  in  uneven  resis¬ 
tance  increase  and  degradation  in  the  cell,  which  may  result  in 
temperature  increase  and  even  thermal  runaway.  Fig.  13  shows  the 
capacity  loss  due  to  SEI  growth.  The  capacity  loss  stabilizes  when 
the  SEI  layer  reaches  its  stable  thickness.  The  SEI  layer  reaches  its 
stable  state  fairly  quickly  due  to  the  slowdown  of  the  diffusion  and 
kinetic  reaction  processes,  which  is  consistent  with  several  ob¬ 
servations  [15,40,43  . 

4.  Conclusions 

In  this  study  a  cell-level  physics-based  model  has  been  devel¬ 
oped  to  couple  the  SEI  layer  growth  and  the  thermal  effect.  The 
model  provides  an  effective  tool  to  investigate  and  evaluate  the 
effects  of  SEI  layer  growth  and  temperature  increase  on  the  per¬ 
formance  of  LIBs.  Depending  on  the  processes  of  diffusion  and  side 
reaction,  the  SEI  layer  grows  quickly  to  the  stable  state  in  the  first 
several  cycles.  The  slowing  down  of  the  SEI  layer  growth  can  be 
attributed  to  the  difficulty  in  the  diffusion  process.  The  SEI  growth 
under  charge  and  discharge  conditions  shows  different  rates  due  to 
the  electron  flux  difference  through  the  SEI  layer  and  temperature 
change  during  cycling.  The  SEI  growth  rate  shows  position 
dependence  in  the  cell  since  the  current  flow  and  lithium  diffusion 
vary  at  different  locations  in  the  electrode.  The  SEI  grows  faster  at 
the  interface  between  the  anode  and  the  separator,  which  is  the 
strategic  location  for  battery  health  monitoring  in  the  battery 
control  algorithm.  The  effect  of  environmental  temperature  and  C- 
rate  on  SEI  layer  growth  has  also  been  investigated  in  this  study.  In 
the  future,  we  plan  to  expand  the  model  to  perform  EIS  simulation, 
and  to  couple  additional  degradation  mechanisms,  such  as  lithium 
plating  and  manganese  deposition.  The  model  framework  can  be 
used  as  an  effective  tool  for  designing  and  optimizing  the  cell 
performance. 
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